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Abstract 

We show that various models of the elliptic Calogero-Moser systems are accom- 
panied with an isomonodromic system on a torus. The isomonodromic partner is a 
non-autonomous Hamiltonian system defined by the same Hamiltonian. The role of 
the time variable is played by the modulus of the base torus. A suitably chosen Lax 
pair (with an elliptic spectral parameter) of the elliptic Calogero-Moser system turns 
out to give a Lax representation of the non-autonomous system as well. This Lax 
representation ensures that the non-autonomous system describes isomonodromic 
deformations of a linear ordinary differential equation on the torus on which the 
spectral parameter of the Lax pair is defined. A particularly interesting example 
is the "extended twisted BCi model" recently introduced along with some other 
models by Bordner and Sasaki, who remarked that this system is equivalent to In- 
ozemtsev's generalized elliptic Calogero-Moser system. We use the "root type" Lax 
pair developed by Bordner et al. to formulate the associated isomonodromic system 
on the torus. 
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1 Introduction 



In 1996, Manin U proposed a new expression of the sixth Painleve equation. This is a 
differential equation of the form 



12, 3 



(27u) 2 ^ = 5> aP '(g + ^), (1.1) 



dr 2 

a=0 

where p'(u) is the derivative of the Weierstrass p function with primitive periods 1 and 
r. 



p(u) = p(u 
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u 2 ^— ' \ (u + m + nr) 2 (m + nr) 2 J ' 

(m,n)^(0,0) VV 1 V ' 7 

u; a (a = 0, 1, 2, 3) are the origin and the three half-periods of the torus E T = C/(Z + rZ), 

1 It r 

w = 0, wi = -, w 2 = - + -, ^3 = ^' ( L3 ) 

and a a (a = 0, 1, 2, 3) are the simple linear combinations (aco, «i, «2, «3) = (a, — /3, 7, 1/2 — 
5) of the four parameters a, /3, 7 and /3 of the sixth Painleve equation 

^! = 1 p 1 1 1 1 ) ( d y) 2 pi 1 1 1 ) d v 

dx 2 2\y y — 1 y — x J \dx J \x x — 1 y — x J dx 

■u(-u — l)(w — x) / a; — 1 r x(x — 1)\ .„ 

Manin's equation can be written in the Hamiltonian form 

dq dp &H 

2m Tr =P > 2m Tr = ~^ (L5) 

with the Hamiltonian 

1 3 

H = -p 2 - ^2a a p(q + u a ). (1.6) 

a=0 

Since the Hamiltonian depends on the modulus r explicitly, this is a non-autonomous 
Hamiltonian system. In this new framework, Manin reconsidered the affine Weyl group 
symmetries of the sixth Painleve equation discovered by Okamoto , solutions for special 
values of a, /3, 7 and 5 constructed by Hitchin ||, etc. 

Manin's equation reveals an unexpected link between the Painleve equation and the 
elliptic Calogero-Moser systems, i.e., the Calogero-Moser systems 0] with elliptic poten- 
tials. In order to see this relation, we introduce a new variable t and formally replace 
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2nid/dr — > d/dt in the aforementioned equations. The outcome are the autonomous 
equation 

d 2 q 3 

= ^2a a p'(q + u a ) (1.7) 

a=0 

and its Hamiltonian form 

dq dp dH 

Tt =p ' Tt = ~W (L8) 

If all a n 's take the same value —g 2 /8, one can use an identity of the p function to rewrite 
the above equation as: 

^ = -ji2p'(<l + ^) = -9 2 p'(^)- (1-9) 

a=0 

This is exactly the two-body elliptic Calogero-Moser system; the £-body elliptic Calogero- 
Moser system (v4^_i model) is defined by the Hamiltonian 

1 £ 2 

w = 2&?+ 7 !>(«»■-»)■ ( L1 °) 

j=l j^k 

As Krichever || demonstrated, this elliptic Calogero-Moser system is an isospectral inte- 
grable system with a Lax representation 

^- = [L(z),M(z)\, (1.11) 

where the Lax pair L(z) and M(z) are matrix-valued functions of a spectral parameter z 
on the torus E T . Furthermore, the general case falls into Inozemtsev's generalization of 
the elliptic Calogero-Moser system || defined by the Hamiltonian 

n = \ s trt + q f E E^+ e ^) + ^EE^fe+^)- (i.i2) 

3=1 e,e'=±l j^k j=l a=0 

Levin and Olshanetsky [fjj developed a geometric formulation of isomonodromic sys- 
tems on a general Riemann surface, and characterized Manin's equation as an isomon- 
odromic system on the torus E T . Their interpretation of isomonodromic deformations 
is based on the notion of the Hitchin systems |fl. According to this interpretation, the 
coordinates qj of Calogero-Moser particles are identified with the moduli of an SU(£) flat 



bundle on the torus E T , and the L-matrix L(z) is nothing but the Higgs field on this bun- 
dle. (Such a link between the elliptic Calogero-Moser systems and the Hitchin systems 
was already pointed out before their work by Nekrasov [|9] and Enriquez and Rubtsov 
1(J.) Isomonodromic deformations are special deformations of these geometric data as 



the complex structure of the base torus (or, equivalently, the modulus r) varies. This ge- 
ometric picture suggests a wide range of generalizations of isomonodromic deformations 
(see, e.g., the recent work of Levin and Olshanetsky ||1 1|| ) . 

Unfortunately, however, it is only the special case with ao = cti = a 2 — ^3 that was 
successfully treated in the formulation of Levin and Olshanetsky. This is simply because 
no suitable Lax representation was available for the Inozemtsev system. Inozemtsev |J 
presented a Lax representation, but it is not suited for that purpose. 

Recently, a new type of Lax pair — the root type Lax pair — was proposed by Bordner 
et al. fl2| , |13| , |14f for various models of the elliptic Calogero-Moser systems including the 
Inozemtsev system. This is a Lax pair constructed on the basis of an underlying root 
system (e.g., the Ag^\ root system for the aforementioned elliptic Calogero-Moser system, 
and the BCi root system for the Inozemtsev system). The construction covers not only the 
ordinary elliptic Calogero-Moser systems (the "untwisted models" ) but also the "twisted 
models" introduced by D'Hoker and Phong |15|] and their generalizations (the "extended 
twisted models"). The Inozemtsev system coincides with the extended twisted BCe model 
in the classification of Bordner and Sasaki []T4" |. In particular, the root type Lax pair for the 
extended twisted BC\ model gives a Lax representation to the aforementioned isospectral 
analogue of Manin's equation. 

One of the goals of this paper is to show, using the root type Lax pair, that each of 
these elliptic Calogero-Moser systems are accompanied with an isomonodromic system on 
a torus. The fist step of the construction is simply to replace the equations of motions 

* = {,,„}, * = feH} (1 ,3) 

of the elliptic Calogero-Moser system by the non-autonomous system 

2 m ^ = {q,H}, 2m^ = {p,H} (1.14) 

with the same Hamiltonian TC. We then rewrite this non-autonomous system into a Lax 
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equation of the form 



M ?m + 8w = miU(z)] (1 . 15) 



using a root type Lax pair L(z) and M(z). This Lax equation implies the Frobenius 
integrability of the linear system 

^l = L(z)Y(z), 2 m ^- + M(z)Y(z)=0, (1.16) 

from which one can deduce that the non-autonomous system is an isomonodromic system 
on the torus E T . 

Actually, we shall use the root type Lax pair made of slightly different building blocks. 
The root type Lax pairs, like the previously known Lax pairs, contain complex ana- 
lytic functions x(u,z), y(u,z), etc. that satisfy special functional equations (called the 



"Calogero functional equations" ||16|| ). Bordner et al. use the Weierstrass sigma function 
to construct those functions. We use the Jacobi theta function 9\ instead. This is inspired 
by the work of Levin and Olshanetsky, who used substantially the same function to con- 
struct the L- matrix (i.e., the Higgs field in their framework) for isomonodromic systems 
on a torus. This minuscule difference is rather crucial for deriving an isomonodromic Lax 
equation as above. 

The functions x(u, z) and y(u, z) that we use are, in fact, identical to the functions 
that Felder and Wieczerkowski [|17| used in their study on the Knizhnik-Zamolodchikov- 
Bernard (KZB) equation This is by no means a coincidence. As Levin and Olshanet- 
sky stressed, the KZB equation and the Hitchin system (or, rather, its isomonodromic 
version) are closely related. 

In order to illustrate that our method also works for some other cases, we show a 
construction of an isomonodromic analogue for the "spin generalization" [|19| of the elliptic 
Calogero-Moser system. Actually, a multi-spin generalization of this construction is also 
possible, which is nothing but the genus-one case of Levin and Olshanetsky's framework. 

This paper is organized as follows. In Section 2, we illustrate our construction of 
isomonodromic systems in the case of the most classical ^Lj_i model. This will serve as a 
prototype of the subsequent discussion. Section 3 is devoted to the models treated by the 
root type Lax pairs, and Section 4 to the spin generalization. Section 5 is for concluding 
remarks. Technically complicated calculations are collected in Appendices. 
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2 Isomonodromic Systems on the Torus — a Proto- 
type 

We start with illustrating our construction for the most fundamental case — the the A^_i 
model and its Lax pair in the vector representation of SU(£). 

2.1 Ai-i Model of Elliptic Calogero-Moser Systems 

The Ag_i model is defined by the Hamiltonian 

3=1 j^k 

Here qj and pj (j = 1, • • • ,£) are the coordinates and momenta of the particles with the 
canonical Poisson brackets 

{<lj,Pk} = 5jk, {qj, q k } = {Pj,Pk} = 0. (2.2) 

Following Manin's equation, we noralize the primitive periods as 

2cui = 1, 2lu 3 = r (2.3) 

The equations of motion are give by the canonical equations 

dqj , , 

~dt = te' W >=*' 

^ = {Pj,H} = -J , Ep'(«'-ft)- ( 2 - 4 ) 
This elliptic Calogero-Moser system has a Lax pair of the form 

L ( z ) = ^PjEjj + i 9^2 x ^ ~ qk > z ) E i^ 

3=1 j¥=k 
I 

M(z) = Y, D j E 3j J rwY.y^-^z)E jk , (2.5) 

3=1 j¥=k 

where Ej k is the matrix unit, (Ej k ) mn = S m j5 nk . The diagonal elements Dj of M(z) are 
given by 

D j = i 952 fp ( q i~ qk ')> ( 2 - 6 ^ 
6 



and x(u, z) is a function that satisfies, along with its w-derivative 



the functional equations 



x(u, z)y(v, z) - y(u, z)x(v, z) = x(u + v, z)(p(u) - p(v)), (2.8) 
z)y(-u, z) - y(u, z)x(-u, z) = p'(u), (2.9) 
x(u, z)x(—u, z) — p(z) — p(u). (2-10) 

Using these functional equations, one can easily prove the following well known result || : 

Proposition 1 The matrices L(z) and M(z) satisfy the Lax equation 

^=[L(z),M(z)}. (2.11) 

As far as the elliptic Calogero-Moser system is concerned, the choice of x(u, z) and y(u, y) 
is rather irrelevant. A standard choice is the function 

o~(z — u) , 
x(u,z)= , 2.12 

a{z)a{u) 

where <j(u) = a{u \ l,r) is the Weierstrass sigma function with primitive periods 1 and 

T. 

Thus, the elliptic Calogero-Moser system is an isospectral integrable system. An 
involutive set of conserved quantities can be extracted from the traces Tr L(z) k , k = 
2, 3, • • • of powers of the L-matrix. The quadratic trace is substantially the Hamiltonian 
itself: 

L(z) 2 

Tr — - — = H + (independent of p and q). (2.13) 

The functions x(u, z) and y(u, z) based on the sigma function, however, are not very 
suited for constructing an isomonodromic system. We shall show an alternative in the 
next subsection. 
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2.2 Our choice of x(u, z) and y(u, z) 

Inspired by the work of Levin and Olshanetsky J7|, we take the following function x(u, z) 
and its w-derivative y(u, z) for constructing an isomonodromic Lax pair: 

xM ^j^m. (2 , 4) 



exp Lit (n+^j + 2vrz (n+ (u+ V (2.15) 



Here 9\(u) is one of Jacobi's elliptic theta functions, 
9 1 (u) = 9 1 (u\t)-- 

and 9[(u) its derivative. Accordingly, the partner y(u,z) can be written 

y(u, z) = -x(u, z)(p(u) + p(z - u)), (2.16) 
where p(u) denotes the logarithmic derivative of 9\(u), 

Pin) = (2.17) 
The function p(w), too, plays an important role throughout this paper. 



Proposition 2 These functions x(u, z) and y(u, z) satisfy the functional equations {\2.8j ) 
- ([2.10] ) and the differential equation 

2 m d -^ + ^^l = 0. (2.18) 
or ouaz 

The last differential equation (a kind of 1 + 2-dimensional "heat equation" ) is a char- 
acteristic of our (x, y) pair, and plays a key role in our construction of isomonodromic 
systems. 

We give a proof of these properties in Appendix A. The following are supplementary 
remarks on these functions. 



The proof of ( [2.8H2.10D is based on the following analytical properties of x(u, z): 



1. x(u,z) is a meromorphic function of u and z. The poles on the u plane and 
the z plane are both located at the lattice points u = m + nr and z = m + nr 
(m, n £ Z). 
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2. x(u, z) has the following quasi-periodicity: 

x{u + 1, z) = x(u, z), x{u + r, z) = e 2mz x{u, z), 

x(u, 2 + 1) =x(u, z), x(u,z + t) = e 27riu x(u, z). (2.19) 

3. At the origin of the u and z planes, x(u, z) exhibits the following singular 
behavior: 

x(u,z) = p(z)+0(u) (u—>0), 

u 

x(u,z) = -- + p{u) + 0{z) (z->0). (2.20) 

z 

These properties are an immediate consequence of the following well known fact: 

1. 6\(u) is an entire function with simple zeros at the lattice points u = m + nr 
(m, n G Z). 

2. 6i{u) is an odd and quasi-periodic function, 

9 1 (-u) = e 1 (u + i) = -e 1 (u), 

9 x {u + r) = -e- mT - 2mu Q x [u). (2.21) 

One can similarly see the following analytical properties of p(u): 

1. p(u) is a meromorphic function with poles at the lattice points u = m + nr 
(m, n G Z). 

2. p(w) is an odd function with additive quasi-periodicity: 

p(-u) = -p{u), p(u + l) = p{u), p{u + r) = p(u) - 2m. (2.22) 

3. At the origin u = 0, p(tt) exhibits the following singular behavior: 

P^) = - + ^-f + 0(u 3 ) (u-.0). (2.23) 



The proof of ( [2.18 ) is based on the well known "heat equation" 



= ^iM. (2.24) 



of the Jacobi theta function. 



2.3 Isomonodromic deformations 



Replacing d/dt — ► 2md/dT, one obtains a non- autonomous Hamiltonian system: 

2™^ = fe, H} = -g 2 p'fa ~ ft)- (2-25) 

We now demonstrate that this gives an isomonodromic system on the torus E T . A key is 
the following Lax equation: 

Proposition 3 L(z) and M(z) satisfy the Lax equation 

J^+^ = [K,),«(4 (2.26) 

Proof. Let us notice that the right hand side of the isospectral Lax equation is in fact the 
Poisson bracket of L(z) and the Hamiltonian: 

[L(z),M(z)] = ^ = {L(z),H}. (2.27) 

Since the phase space and the Hamiltonian are the same as those of the original system, 
the relation [L(z), M(z)] = {L(z),TC} persists in the present setup. Thus the right hand 
side of the Lax equation can be written 

[L{z),M(z)) = {L(z),H} 
t 

= X^'' H } E n + W ^2{Qj - Qk, T-L}y(qj ~ Qk, z)E jk . (2.28) 

j=l j^k 

On the other hand, 

dL(z) dM(z) d Pj „ 



+i g J2 2ni {^~ ~d~) yfa - qk ' z "> E J k 

+!g W 2 „^+*^£T) Ej , (2 . 29) 

j^k ^ / u=q 3 -q k 

The last sum vanishes because of the "heat equation" (|2.18|) . The other part coincides, 
term-by-term, with the above expression of the commutator [L(z), M(z)]. Q.E.D. 
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This Lax equation enables us to interpret the non-autonomous Hamiltonian system as 
an isomonodromic system on the torus E T . The Lax equation is nothing but the Frobenius 
integrability condition of a linear system of the form 

= L(z)Y(z), 2™^^ + M(z)Y(z) = 0. (2.30) 

The first equation is an ordinary differential equation on the torus E T , and has a regular 
singular point at z — 0. Analytic continuation of the solution around this singular point 
yields a monodromy matrix To- Besides this local monodromy matrix, there are global 
monodromy matrices T a and Tp that arise in analytic continuation along the a (z — > z+1) 
and (3 (z — > z + r) cycles. The second equation of the above linear system implies that 
these monodromy matrices are left invariant as r varies. 

Let us specify this observation in more detail. The situation is more complicated 
than isomonodromic systems on the Riemann sphere: The monodromy of L(z) and M(z) 
themselves are non-trivial, 

L(z + l)=L(z), M(z + l)=M(z), 
L(z + r) = e 2mQ L(z)e- 2nlQ , 

M(z + r) = e 2niQ (M(z) + 2mL(z))e^ iQ - 2niP, (2.31) 

where Q = Y2j=iQjEjj and P = Ylj=iPj-^jj- These relations are a consequence of the 
quasi-periodicity of x(u,z), y(u,z) and p(z). The monodromy of L(z) implies that Y(z) 
has to be treated as a section of a non-trivial GL(£, C)-bundle (or SL (£, C)-bundle, if 
we take the center of mass frame with Yl]=iPj = 0) 011 ^ ne torus E T . The monodromy 
matrices r , T a and Tp thus arise as follows: 

Y (ze 2m ) = Y(z)T , Y(z + 1) = Y(z)T a , Y(z + r) = e 2mQ Y{z)T p . (2.32) 

Note that the exponential factor in the last relation reflects the non-trivial monodromy 
of L(z) along the /3-cycle. Having this monodromy structure of Y(z), one can deduce the 
following fundamental observation: 



Proposition 4 The monodromy matrices do not depend on r, i.e., 

= 0. (2.33) 



dT _ dT a _ dYp 



dr dr dr 
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Proof. Let us rewrite the second equation of the linear system as 



M(z) = -2ni 



dY(z) 



Y(z) 



-i 



(2.34) 



and examine the implication of the monodromy structure of Y(z) noted above. This leads 
to the following relations: 



(We have used the relation 2-KidQ/dr = P.) These relations are consistent with the 
aforementioned monodromy structure of M(z) if and only if the monodromy matrices of 
Y(z) are independent of r. Q.E.D. 

3 Elliptic Calogero-Moser Systems Based on Root 
Systems 

Here we consider the elliptic Calogero-Moser systems associated with a general irreducible 
(but not necessary reduced) root system A. 

In the following, the root system A is assumed to be realized in an ^-dimensional 
Euclidean space M — M. . Let x ■ y denote the inner product of two vectors in M and 
its bilinear extension to the complexification M c = M £g> K C. The dual space M* = 
Hom(M, R) of M is identified with M by this inner product. Each element a E A induces 
a reflection (the Weyl reflection) s a (x) = x — (2a ■ x/a ■ a)a. This gives a representation 
of the Weyl group W(A) on M. The root system A is invariant under the action of this 
Weyl group. 

The elliptic Calogero-Moser system associated with the root system A is a Hamiltonian 
system on M x M (or its complexification M c x M c ). The orthognal coordinates (q,p) = 
(qi, • • • , qe,Pi, • • • ,pe) of M x M give canonical coordinates and momenta with the Poisson 



M(z+1) 
M(z + r) 



M(ze 2ni ) 




(2.35) 
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brackets 



{Qj,Pk} = 8jk, {Pj,Pk} = {qj,qk} = 0. (3.1) 
3.1 Simply laced models 

We first consider the case of simply laced (Ae-i, De and Eg) root systems. The associated 
elliptic Calogero-Moser system is defined by the Hamiltonian 

Here g is a coupling constant, and p(u) the Weierstrass p function with primitive periods 
1 and r. The equations of motion can be written 

^=P, ^ = -j^p(a. ff )a. (3.3) 

aSA 



We first review the "root type" Lax pair of Bordner et al. for these models |T|], then 
explain how to convert these isospectral systems to isomonodromic systems. 

3.1.1 Root type Lax pair 

The "root type" Lax pair for these simply laced models are A x A matrices, i.e., matrices 
whose rows and columns are indexed by the root system A. They are made of three parts: 

L(z) =P + X l {z) +X 2 (z), M{z) = D + Y 1 {z) + Y 2 {z). (3.4) 
P and D are diagonal matrices, 

Pp~f=P- /3<W D Pl = ZV/37 (P> 7 e A), (3.5) 
and the diagonal elements Dp of D are given by 

Dp = igp((3-q)+ig p(l ■ q)- (3.6) 

7GA,/3-7=l 

Xi(z), etc. are diagonal- free matrices of the form 

X\{z) = ig x(a ■ q, z)E(a), 

oSA 
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X 2 (z) = 2ig^2x{a-q,2z)E(2a), 
Y i( z ) = ig^2v(,a-q,z)E(a), 

Y 2 {z) = ig^2y(a-q > 2z)E(2a) > (3.7) 

aGA 

where x(u, z) and y(u, z) are the same as the functions used in the previous section, and 
E(a) and E(2a) are A x A matrices of the form 

E(ai)pj = 5 Q ,/3- 7 , E(2a)p 7 = 5 2 a,/?- 7 (A7GA). (3.8) 

(We have slightly modified the notation of Bordner et al: x(u,2z), y(u,2z) and E(2a) 
amount to Xd(u, z), yd(u,z) and Ed(a) in their notation.) 
These matrices satisfy the Lax equation 

^ = [L(z),M(z)} (3.9) 

under the equations of motions. The traces Tr L(z) k , k = 2,3,---, of powers of L(z) 
are conserved, and an involutive set of conserved quantities can be extracted from these 
traces. The Hamiltonian itself can be reproduced from the quadratic trace TrL(z) 2 . We 
refer the details of these results to the paper of Bordner et al. [12|]. The choice of x(u, z) 
and y(u, z) is irrelevant in this case, too. 

Thus, in particular, the A^_i model turns out to have at least two distinct Lax pairs 
- the Lax pair of £ x I matrices realized in the vector representation of sl(£), and the 
Lax pair of £(£ — 1) x £(£ — 1) matrices based on the A^i root system. This is also the 
case for the other simply laced root systems. Bordner et al. call the Lax pairs of the first 
type the "minimal type", because they are realized in a minimal representation of the 
associated (not necessary simply laced) Lie algebra. It should be noted that the "root 
type" Lax pairs do not possess a Lie algebraic structure; unlike the usual root basis of 
simple Lie algebras, the matrices E(a) and E(2a) are not closed under the Lie bracket. 

3.1.2 Isomonodromic system 

The prescription for constructing an isomonodromic analogue is the same as the previous 
case, namely, to replace d/dt — > 2irid/dT. This converts the equations of motion of the 
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elliptic Calogero-Moser system to the non-autonomous system 



aeA 



Let x(u, z) be the function defined in Q2.14Q , and y(u,z) its w-derivative. The following 



are the keys to an isomonodromic interpretation. 



Proposition 5 1. L(z) and K(z) satisfy the Lax equation 



2 J_m + amA = [L(z)>M{z)] (311) 

2. L(z) and M(z) have the following monodromy property: 

L(z + l)=L(z), M(z + l) = M(z), 
L(z + t) = e 2niQ L(z)e- 2mQ , 

M{z + t) = e 2wiQ {M{z) + 2mL{z))e- 2wiQ - 2mP, (3.12) 

where Q is the diagonal matrix with matrix elements = q ■ (38 p^. 

Proof. The proof is almost the same as the proof for the isomonodromic Lax pair of the 
Ai-i model in the vector representation. Let us first verify the Lax equation. The right 
hand side of the Lax equation can be written 



[L(z),M(z)\ = {P,H} + ig 1 £ t {a-q,n}y(a-q,z)E( a ) 

aeA 

+2ig ^{a ■ q, H}y{a ■ q, 2z)E(2a). (3.13) 



On the other hand, 



„ .dL(z) dM(z) n .dP . ^ n .da-q , , _ , 

aeA 

+2ig 2m-^-y(a ■ q, 2z)E(2a) 

aeA 

+ls W 2m ^ + M-£>) E{a) 

aeA ^ ' u=a-q 

(, .dx(u,2z) dy(u,2z)\ . . „ 

+2ig 4 ^^H + E{2a) - (3 ' 14) 

aeA ^ ' u=a-q 
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The last two sums vanish because of Q2.18| ). The other part coincides by the equations of 
motion. Thus we obtain the Lax equation. Let us next consider the monodromy of L(z) 
and M(z). Note the commutation relations 

[Q,E(a)]=q-aE(a), [Q,E(2a)) = 2q-aE(2a), (3.15) 

which can be exponentiated as follows: 

e 2wiQ E(a)e- 2mQ = e 2mq - a E(a), e 2mQ E(2a)e~ 2mQ = e 4mg - a E(2a). (3.16) 

The monodromy property of L(z) and M(z) can be derived from these relations and the 
quasi-periodicity of x(u, z) and y(u, z). Q.E.D. 

The rest is parallel to the case in the previous section. The only difference is that the 
ordinary differential equation 

d ^ = L{z)Y{z) (3.17) 

on the torus E T has four regular singular points at z — 0, co>2, uj^. The latter three 
singular points originates in Xz(z). Let r a (a = 0, 1, 2, 3) denote the monodromy matrices 
in analytic continuation of Y(z) around these four points. The Lax equation implies that 
these local monodromy matrices and the two global ones T a and Y p are independent of r: 

0. (3.18) 



Or dr Or Or 



3.2 Non-simply laced models 

The elliptic Calogero-Moser system associated with a non-simply laced (B^, Ce, G<i and 
BCe) root systems can have several independent coupling constants, one for each Weyl 
group orbit in the root system. The root type Lax pairs are extended to the non-simply 
laced cases by Bordner et al. JTBJ. As they pointed out, one can construct a different root 
type Lax pair for each Weyl group orbit of the root system. Thus the B £ , Ci, F 4 and G 2 
models have, respectively, two distinct Lax pairs based on the orbits of long and short 
roots, whereas the BCi model has three based on the orbits of long, middle, and short 
roots. Note that each Weyl group orbit consists of roots of the same length. 

Although all the non-simply laced models can be treated in the same way, let us 
illustrate our construction of isomonodromic systems for the BCi model. This is also 
intended to be a prototype of the case that we shall consider in the next subsection. 
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3.2.1 BC e model 

The BCg ro °t system can be realized in M = M 1 : 



A(BC e ) = A ; UA m UA S) 

A t = {±2ej | 1 < j < £} (long roots), 

A m = {±ej ± et | j 7^ k} (middle roots), 

A s = {±ej | 1 < j < £} (short roots), (3.19) 

where ei, • • • , eg are the standard orthonormal basis of ~R e . A;, A rn and A s give the three 
Weyl group orbits. 

The Hamiltonian of the BCe model takes the form 

n = \w+ 9 \ ?(«•?) + j£p(«-9)+& 2 Z>(«-9)- (3- 2 °) 

a€A m aeA; aeA s 

The equations of motion can be written 
dq 

TT = V \ 

% = - g -f E ■ 9)« - f E • 9)« - ^ 2 E • ^)«- ( 3 - 21 ) 

aeA m aeA ; aeA s 

g m , gi and g s are three independent coupling constants. g s is a modified ( "renormalized" in 
the terminology of Bordner et al.) coupling constant connected with a more fundamental 
( "bare" , so to speak) coupling constant g s as 

~9 2 s = 9 2 s + 9 ~f- (3-22) 
The "bare" coupling constant appears in the construction of a Lax pair. 



3.2.2 Root type Lax pair for BCi model 

As mentioned above, there are at least three root type Lax pairs based on the three Weyl 
group orbits A m , A^ and A s . Bordner et al. constructed only one of them, namely, a Lax 
pair based on A m . Here we present a Lax pair based on A s . This is a 21 x 2£ system, 
much smaller than the Lax pair based on A m , and presumably more suitable for studying 
the associated isomonodromic deformations. 
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The Lax pair are indexed by A s and take the following form: 

L(z) = P + X 1 (z)+X 2 (z)+X 3 (z), 
M(z) = D + Y 1 (z)+Y 2 (z)+Y 3 (z). 

P and D are diagonal matrices, 

Pp~<=v- AW D Pi = D P s /3~/ (P, 7 e A s ), 
and the diagonal elements of .D are given by 

Dp = ig m p(l- ( l)+ i 9ip(^P-(l)+Wsp((3-q). 

7 eA s ,/3-7=l 

Xl(z), etc. are diagonal- free matrices of the form 

X 1 (z) 



X 2 (z 
X 3 (z 
Y 1 (z 
Y 2 (z 
Y 3 (z 



= Wm ^ X ( a ' Q> Z ) E ( a )> 

= igi^2 i x(a- q,z)E(a), 

aeAi 

= 2ig s x(a-q,2z)E(2a), 

aeA s 

= Wm ^ y ( a ' Q> z ) E ( a )> 
= igi^2y(a-q,z)E(a), 

= Ws ^2 y(a ■ q,2z)E(2a), 



aeA s 



where 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



(3.27) 



E(a) Pl = <J Q>/ 3_ 7 , E(2a)fr = 5 2 a,/3~ 7 (/?, 7 e A s ). 

This Lax pair is a specialization of the Lax pair for the extended twisted model that we 
shall present in the next subsection. 



3.2.3 Isomonodromic system 

This system, too, can be converted to an isomonodromic system by replacing d/dt 
2nid/dr. The equations of motion are a non-autonomous system of the form 

2m— = p, 
dr 



2ni^- = - 9 -f Y, P'(<* ■ 9)« - f E ^ a ■ ^ a - % E • (3-28) 



2 ^— ' " 4 

QGA m oeA ; aeA s 
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The following can be verified just as in the case of simply lased models: 

1. L(z) and M(z) satisfy the Lax equation 

2«^ + ^M = [L W ,M W ]. (3.29) 

2. L(z) and M(z) have the following monodromy property: 

L(z + l) = L(z), M(z + l) = M(z), 
L(z + t)= e 2niQ L(z)e- 2mQ , 

M(z + r) = e 2niQ (M(z) + 2niL(z))e- 2niQ - 2niP. (3.30) 

The interpretation of this Lax equation, too, is parallel to the simply laced models. The 
ordinary differential equation 

^ = L{z)Y{z) (3.31) 
dz 

on the torus E T has four regular singular points at z — 0, Ui, u>2, 003. The local monodromy 
matrices T a (a = 0, 1, 2, 3) at these points and the global monodromy matrices T a and Tp 
are invariant as r varies. 

3.3 Twisted and extended twisted models 

We now proceed to the "twisted" and "extended twisted" models. The Hamiltonian of 
the untwisted models can be generally written 

K= l -V-V+\Y.9 2 \ a \^ a - ( l)- ( 3 - 32 ) 

aeA 

The twisted models, introduced by D'Hoker and Phong \T5\ for non-simply laced root 
systems, are defined by a Hamiltonian of the form 



'H^-p-p+-^2gf a] p u ( a )(a-q), (3.33) 

where p u (a){ u ) are the p-functions with suitably rescaled primitive periods. D'Hoker 
and Phong proved the integrability of those twisted models by constructing a Lax pair 
in a representation of the associated Lie algebra. Bordner and Sasaki |14| proposed an 
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alternative approach based on root systems rather than Lie algebras, and pointed out that 
the twisted model of the B#, Ce and BCi types can be further extended. The extended 
twisted models have one (for the Be and Ce models) or two (for the BCe model) extra 
types of elliptic potentials. 

Our construction of isomonodromic systems can be extended to the twisted and ex- 
tended twisted models. We illustrate this result, just as in the previous subsection, for 
the BCe model. As Bordner and Sasaki noted, the extended twisted BCe model is made 
of five different types of elliptic potentials, and coincides with the Inozemtsev system ||. 

3.3.1 Extended twisted BCi model 

The extended twisted BCe model is defined by the Hamiltonian 

aeAm aeAi aeA; 

+fsi E p(« • 9) + & E ^ (1/2) (« • «)■ ( 3 - 34 ) 

a£A s agA s 

§12, 9 si and g s2 are "renormalized" coupling constants, which are related to unrenormalized 
coupling constants gi 2 , g s i and g s2 as follows: 

9i2 = 9h + 2 9n9i2, 

9 si = 9si + 2 9si9s2 + -jiOsigix + g s i9i2 + 9,2912) , 

~2 2 , 9s29ll /„ QC .x 

9s2 = &2 + — g— • ( 3 - 35 ) 

p^ 1 / 2 ) and are the p functions with rescaled primitive periods: 

p^ 2 \u) = p(u I It), p (2 \u) = p(u | 2,r). (3.36) 

(This Hamiltonian is slightly different from the Hamiltonian of Bordner and Sasaki, 
though the contents are essentially the same. With this modification, this model reduces 
to the untwisted BCe model as g\ 2 — * and g s2 — ► 0.) 

3.3.2 Root type Lax pair for extended twisted BCe model 

One can construct, like the untwisted model, three different root type Lax pairs can be 
constructed based on the three Weyl group orbits A m , A/ and A s . The Lax pair based 
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on A m is presented by Bordner and Sasaki. The Lax pair based on A s can be obtained 
by modifying the Lax pair for the untwisted BCi model as follows. 

The Lax pair L(z) and M(z) are indexed by A s and made of four parts, 

L(z) = P + X 1 (z)+X 2 (z)+X 3 (z), 

M(z) = D + Y 1 {z)+Y 2 {z)+Y 3 {z). (3.37) 

The diagonal matrix P is the same as the P in the untwisted model. The diagonal 
matrices of D are given by 

Dp = ig m p ^ ' q "> + ^iipC 2 ^ • q) + i 9i2p {2) (2(3 • q) 

7eA m ,/3- 7 =l 

+ig sl p(f3-q)+tg s2 p^(f3-q). (3.38) 

Xi(z) and Yi(z) are the same as those for the untwisted model. The other matrices take 
the following form: 



X 2 {z) = ^ (igixx{a -q,z)+ igi 2 x® (a ■ q, z) j E(a) , 



hig s ix{a ■ q, 2z) + 2ig s2 x {1/2) (a ■ q, 2z)\E{2a), 



oSA 3 

Y 2 (z) = ^2^igi 1 y(a-q,z)+ig l2 y (2) (a-q,z)^E(a), 

Ys(z) = ^2(ig sl y{a-q,2z)+ig s2 y^ 2 \a-q,2z))E(2a). (3.39) 

oSA s 

This Lax pair reduces to the Lax pair of the untwisted model if g\ 2 = and g s2 = 0. 

The new objects arising here are the functions x^^ 2 \u, z), x^ 2 \u,z) and their u- 
derivatives 

= »«(«,*) = ^2*. (3.40) 



For the consistency of the Lax equation 

dLjz) 
dt ~ 



[L(z),M(z)}, (3.41) 



these functions have to satisfy several functional equations. D'Hoker and Phong [15] and 



Bordner and Sasaki JTJ] use a set of functions based on the Weierstrass sigma functions. 
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We use the function x(u, z) = x(u, z | r) defined in (|2.14j) and its modifications 

x {1/2) (u,z) = 2x(2u,z\2r) 



20! (z - 


2u 


2t)0' 1 (O 


2t) 


Oi(z\ 


2r] 


\9 x (2u 


2r) 



Oi{z-l\ T -)e'^\ T -) 



2 2 2 2^(z | | I) 



These functions x^ 2 \u, z) and z), too, satisfy 1 + 2-dimensional "heat equations" 

of the form 

dxW>(u,z) d 2 x^(u,z) 
2m a 1 = °' 

OT OUOZ 

2J Xm } U > z) + *^ Z) = 0. (3.43) 

OT OUOZ 

The functional identities for these functions and the proof of the Lax equation are pre- 
sented in Appendices B and C. 

3.3.3 Isomonodromic system 

Replacing d/dt — > 2md/dr, we obtain a non-autonomous Hamiltonian system with the 
same Hamiltonian. The isomonodromic interpretation of this non-autonomous system is 
again based on the following two observations: 

1. L(z) and M(z) satisfy the Lax equation 

2^ + ^M = [L W ,M W1 . (3.44) 

2. The monodromy of L(z) and M(z) is the same as the monodromy of the Lax pair 
for the untwisted model: 

L(z + l) = L(z), M(z + l) = M(z), 
L(z + t)= e 2niQ L(z)e- 2mQ , 

M(z + r) = e 2niQ (M(z) + 2mL(z))e~ 2niQ - 2mP. (3.45) 

The ordinary differential equation defined on the torus E T by the matrix L(z) has four 
regular singular points at u = 0, u\, u 2 , ^3. The Lax equation and the monodromy of L(z) 
and M(z) ensure that the local monodromy matrices T a (a = 0, 1, 2, 3) and the global 
monodromy matrices T a and Tp are independent of r. 
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3.3.4 Relation to Inozemtsev system 

The final task is to clarify the relation to the Inozemtsev system. In terms of the orthog- 
onal coordinates qj — q ■ Cj and pj = p- e 3 - (j — 1, • • • , £), the aforementioned Hamiltonian 
can be written 

j=l e,e'=±l j^k j=l 



+ 9 f J2 P (2 \2 qj ) + 2& J2 pfa) + 2f s2 J2 P (1/2) (^). (3.46) 
j=i j=i j=i 

One can rewrite this Hamiltonian using the identities 

11 1 1 

p{ 2u ) = ^P(u) + -piu + uji) + -p(u + u 2 ) + -p(u + u 3 ), 

P (1/2) {u) = p(u) + p{u + ui) - pM, 

p^\2u) = -^(u)+ l l p{u + u z )-- A p{uz). (3.47) 
The outcome is, up to a term h{r) depending on r only, the Inozemtsev Hamiltonian 

i ,2 ^3 



n = li2p'+ g f E E^+ e ^) + EE^+^) + ^ r )- ( 3 - 48 ) 

j=l e,e'=±l j+k j=l a=0 

The coupling constants g a (a = 0, 1, 2, 3) are given by 



g 2 o = l(g?i + gl) + 2& 2 i + g 2 s2 ), gl = 9 -f + 2^ 2 , 



8 

-.2 



gl = 9 f, gl = l(gl + ~g? 2 )- (3.49) 



4 Spin Generalization of Elliptic Calogero-Moser Sys- 
tems 

"Spin generalization" is a generalization of the elliptic Calogero-Moser systems coupled 
to spin degrees of freedom. Such a spin generalization is characterized by a simple Lie 
algebra rather than a root system. The (classical) spin variables take values in the dual 
space g*, or a coadjoint orbit therein, of the Lie algebra 0. We shall first examine the sl(£) 
model as a prototype, then proceed to the models based on a general simple Lie algebra. 
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4.1 Spin generalization for sl(£) 

The sl{£) spin generalization was first introduced by Krichever et al. |I9| . They obtained 
the spin generalization, just like the spinless case ||, via the pole dynamics of the matrix 
KP hierarchy 

4.1.1 Hamiltonian formalism 

This model is a constrained Hamiltonian system. The Hamiltonian is given by 

1 - 1 

H = 2 ^ ~ 2 ? Pi<lj ~ qk)F ^- 

3=1 j+k 

Here qj and pj (j = 1, • • • , £) are the canonical coordinates and momenta of the Calogero- 
Moser particles, and Fj k (j,k = 1, ■■■,£) a set of classical sl{£) spin variables, whose 
Poisson brackets are determined by the Kostant-Kirillov Poisson structure on the dual 
space of sl(£): 

{Fj k , Fmn\ &mkFjn 3jnF m k- (^"^) 

The equations of motion can be written 
dF- 

= ~ Qm)Fj m + ^ Pfa™ ~ qk)F mk 

-p(Qj - Qk)(Fjj - F kk ). (4.3) 

In particular, the diagonal elements Fjj of the spin variables are conserved quantities: 
dFjj/dt = 0. Although the Hamiltonian does not contain the diagonal elements explicitly, 
they do appear in the equations of motion. We now put the constraints 

F 3J = (j = l, ■■-,£). (4.4) 

These constraints ensure the integrability. (Actually, the integrability is retained if the 
constraints are replaced by Fjj = c, j = 1, • ■ ■ , £, where c is a constant.) 
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4.1.2 Lax pair in vector representation 

The Lax pair of the spinless Ag__\ model in the vector representation of sl(£) can be readily 
extended to the spin generalization as follows: 

t 

L ( z ) = Yl v i E H + Yl a ( q i ~ z ) F *i E 3k, 

j=l j^k 

M i z ) = - /] a (Qj - Qk, z)(p(Qj - gfc) + p(z - Qj + w k ))F kj E jk , (4.5) 

where 

Vi(u) ViizjV^u) 



It is these functions that Felder and Wieczerkowski used in the KZB equation [T^J. The 
function p(u) is already familiar to us. The function a(u, z) is also just a disguise of the 
function x(u, z) that we have used in the preceding sections: 

a(u,z) = —x(u,z). (4.7) 

We however dare to retain the notation of Felder and Wieczerkowski so as to stress the 
similarity with their work. In these notations, the aforementioned functional identities of 
x(u, z) and y(u, z) can be rewritten 

a{u, z)a(v, z)(p(v) + p(z - v) - p{u) - p(z - u)) = a{u + v, z)(p(u) - p(v )),(4.8) 
2a{u, z)a{-u, z){p{u) + p{z - u)) = -p'{u), (4.9) 
a(u, z)a(-u, z) = p(z) - p(u). (4.10) 

Using these functional identities, one can derive the Lax equation 

^ = [L(z),M(z)\. (4.11) 

Note that the constraints ( |4.4j ) are always assumed when we consider the Lax equation. 
Thus the spin generalization, too, is an isospectral integrable system. An involutive set of 
conserved quantities obtained from the traces TrL(z) fc , k — 2, 3, , • • •. The Hamiltonian 
itself can be reproduced from the quadratic trace. 

The matrix F = J2j^k^kjEjk , which is the residue of L(z) at z = 0, stays on 
a coadjoint orbit of sl{€) as t varies. The phase space of the spin generalization can be 
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thereby restricted to the direct product of the phase space of Calogero-Moser particles and 
a coadjoint orbit of various dimensions in the dual space of sl{£). The lowest dimensional 
non-trivial coadjoint orbit can be parametrized by 2£ variables a.j, bj (j — 1, • • • , £) as 

F jk = igb jak (j^k), (4.12) 

where g is a constant. These reduced spin degrees of freedom, however, can be eliminated 
by a diagonal gauge transformation of the Lax equations. (This does not mean that aj 
and bj are non-dynamical. The elimination procedure is done by partially solving the 
equations of motion for those variables.) This gauge transformation in turn gives rise to 
non-zero diagonal elements in M(z), and the outcome is nothing but the Lax equation of 
the spinless elliptic Calogero-Moser system with coupling constant g. The spinless system 
is thus embedded in the spin generalization. 

4.1.3 Isomonodromic system 

There is no substantial difference in the construction of an isomonodromic system. The 
equations of motion are given by 

dqj dpj \ ~\ . 

Mi 

2iri ~^T = ^2p(Qj - Qm)F jm -^2p(q m - q k )F mk . (4.13) 

(Terms including Fjj's have been eliminated by the constraints.) The Lax equation, too, 
can be written in the same form 

n .dL(z) dM(z) rr . . 

Behind this Lax equation is the "heat equation" 

OT ouoz 

satisfied by a(u, z). The final piece of the ring is the monodromy of L(z) and M(z): 

L(z + l)=L(z), M(z + l)=M(z), 
L(z + r) = e 2 ™ Q L(z)e- 2mQ , 

M(z + r) = e 2niQ (M(z) + 2niL(z))e- 2niQ - 2niP. (4.16) 
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As opposed to the root type Lax pairs, the ordinary differential equation 

^ = L(z)Y(z) (4.17) 

on the torus E T has only one regular singularity at z — 0. Thus the local monodromy 
matrix Tq and the global monodromy matrices T a and are all that are invariant under 
the deformations. 

4.2 Preliminaries for general simple Lie algebra 

Let q be a (complex) simple Lie algebra of rank £, f) a Cartan subalgebra, and A the 
associated root system. The Cartan subalgebra induces a root space decomposition of q: 

S = f)©00a- (4.18) 

We choose a basis {e a , \ a G A, (j, — 1, • • • , £} of Q as follows: 

1. /i = !,...,£, are an orthonormal basis of () with respect to the Killing form 
B : fj x 1) -> C, i.e., 

£(/V,M = <W ( 4 - 19 ) 

The Killing form induces an isomorphism ()* = iJom(f),C) ~ f), which determines 
an element /i a for each a G f)*. In terms of the basis of f), this map can be written 
explicitly: 

i 

h a = J2 a (K)K> ( 4 - 2 °) 

2. The root subspace Q a is one dimensional. e a is a basis of g a such that 

[e a ,e_ Q ] = h a . (4.21) 
This choice of e a amounts to the normalization 

B(e a ,e- a ) = l. (4.22) 
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The Lie brackets of the basis elements other than [e a , e 



—a. 



.] now takes the form 



a(/v)e Q , 



0. 



(4.23) 



The structure constants N a ^ are anti-symmetric with respect to the indices, and vanish 
if a + f3 £ A. The following general relation among the structure constants will be used 
in the course of the proof of a Lax equation 



Proof. If a — (3, this relation is trivially satisfied, because all the structure constants 
vanish. Let us consider the case where a ^ (5. By the Jacobi identity, we have 

[e a +p, [e- a ,e-p]\ = [[e a+ p, e_ Q ], e-p] + [e_ Q , [e Q+ ^,e_^]]. 

This implies the identity 



Since we have assumed that a ^ (3, h a and hp are linearly independent, so that the two 
coefficients in this linear retion are equal to zero. Q.E.D. 

We can now specify the classical spin variables for a general simple Lie algebra. Those 
spin variables, by definition, are coordinates of the dual space q* = Hom(Q,C). Let F a 
and be the coordinates dual to the above basis e a and h^. In other words, they are 
the coefficients of e a and hn in the linear combination 



Lemma 1 



N-{3 t a+/3 — N_ a _p — N a+ p_ a . 



(4.24) 



N- a ,-ph a+ p — N a+ p_ a hp — N a+ p_ph, 



which, by the relation h a+ p = h a + hp, can be rewritten 



(N_ a ,_p + N a+ p : _p)h a + (N_ a _p - N a+ p,_ a )hp = 0. 




(4.25) 
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that realizes the isomorphism q* ~ q induced by the Killing form. The Kostant-Kirillov 
Poisson structure on g* determine the Poisson brackets of these spin variables, which take 
the same form as the Lie brackets of the Lie algebra basis: 

i 

{F a ,F„ a } = G a = Y,a(K)G^ 

n=i 

{F a , Fp} = N a:f3 F a+l3 (a + P^O), 
{Gfj,, F a } = a(hn)F a , 

{G„G V } = 0. (4.26) 

4.3 Spin generalization for general simple Lie algebra 
4.3.1 Hamiltonian formalism 

The spin generalization based on q, too, is a constrained Hamiltonian system defined on 
f) x f) x q* by the Hamiltonian 



H = \b(p, P ) - P{u{q))F_ a F a (4.27) 

<*eA 

and the constraints 

G, = (/i=l,---,£). (4.28) 

Here q and p are understood to take values in [). B(p, q) and a(q) amount to p ■ p and 
a ■ q in the models based on root systems. Let us use the same "dot notation" for the 
Killing form f) x [) — > C and the pairing f)* x f) — > C. The Hamiltonian then takes a more 
familiar form: 

n = \p-P-\Y,P(.<X' Q)F- a F a (4.29) 

aeA 

The equations of motion can be readily written down in the language of the coordinates 
Qn — Q • hp and momenta p^ = p ■ of Calogero-Moser particles and the spin variables 
F a and G M on g*: 



dt P ^ 
dv 



dt 2 

aeA 



]- ^2 a ■ h tt<p'{® ■ q)F- a F a , 
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= - Yl q) F a-P F P N a,-p ~ p(a ■ q)G a F a , 

/3eA,o-/3eA 

^ = 0. (4.30) 
dt v 1 

In particular, the diagonal elements of the spin variables are conserved quantities. 
One can thereby safely put the aforementioned constraints. 

4.3.2 Lax pair 

The integrability of our spin generalization is ensured by the existence of a Lax pair as 
follows. 

Proposition 6 Let V be any finite dimensional representation of g, and E a and the 

endomorphisms on V that represent e a and h^. Then the endomorphisms 



a€A £t=l 

M(z) = -J2<7(a-q : z)(p(a-q) + p(z - a ■ q))F„ a E a (4.31) 



a€A 

on V satisfy the Lax equation 

dL(z) 



[L{z),M{z)\. (4.32) 



dt 

Proof. Using the equations of motion and the constraints, one can express the t-derivative 
of the L-matrix as 

dL(z) 



dt 



= 1 + 11 + III, (4.33) 



where 



dt " 2 , 

H=l aeA 



" = EE 



da ■ q da(u, z) 



dt du 

aeA n=l 



F E 

— OL J - J OL 



u=a-q 



= - a ■ aa(a ■ q, z)(p(a ■ q) + p(z - a ■ q))F^ a E a , 

aeA 

III = J2a(a-q,z)^^E a 

oS A 

= - Y o-{a- q,z)p{(3 ■ q)F- a -pFpN- a _pE a . 

Q,/3eA,a+/3^0 
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Similarly, the commutator of the Lax pair can be written 

[L{z),M(z)]=IV + V + VI, (4.34) 
where VI stands for terms from the commutator [P, M(z)], 

IV = - ^ a(a - q, z) (p(a • q) + p(z - a ■ q))F_ a [P, E a ) 
= ~ e(ot ■ q,z)(p(a ■ q) + p(z - a ■ q))a ■ pF„ a E a , 

and V + VI are the the other terms grouped into the Cartan part (V) and the off-Cartan 
part (VI), 

V = - a(-a ■ g, z)a(a ■ q, z)(p(a ■q)+p(z-a- q))F^ a F a [E a , E_ a ) 
- ^ cr(-a ■ q, z)a(a ■ q, z)(p(a ■ q) + p(z - a ■ q))F_ a F a H a , 

aeA 

VI = - v(a-q,q)a(P-q,z)(p(P-q) + p(z-f3-q))F„ a F a {E a ,Ep} 

= - Y o{a ■ q,z)o(l3 ■ q,z)(p((3 ■ q) + p(z - (3 ■ q))F_ a F_pN afi E a+ p. 

It is obvious that IV = II. Using ( |4.9| ), we can readily see that V — I. Thus it remains 
to prove that VI = III. This is achieved as follows: 

vi = -- a ( a ' 9> z ) a (P ' 9, z ) (p(P ■ q) + pO - P ■ q) 

a,/3eA,a+/3^0 

-p(z - a ■ q) - p(a ■ qU F- a F-pN at pE a +p 
[symmetrized with respect to a and /3] 
= ~\ (7 (( a + P) ■ Q,z)(p(a- q) - p(/3 ■ q))F- a F- P N aj(3 E a+ i3. 

a,/3eA,a+/3^0 

[(4.8) is used] 

cr{{a + (3)-q, z)p{(3 ■ q)F_ a F_ p N a ^E a+p 

a,/3GA,a+/3^0 

[asymmetrized with respect to a and 0\ 

= Y a ( a ' 9> z )p{P ■ q)F- a -f3FpN a+ p_pE a . 

a,/3GA,Q+/3^0 

[substituting (3 — > —[3 and a — > a + 0\ 
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Finally using the identity N a+ p-p = —N- a -p, cf. (|4.24j ), we find that the last sum is 
equal to III. Q.E.D. 



Note that the above proof persists to be meaningful if E a and are replaced by the 
Lie algebra elements e a and h^. In other words, the Lax equation actually lives in the 
Lie algebra g itself rather than its representations. This resembles the case of the Toda 
systems. 

4.3.3 Isomonodromic System 

The passage to an isomonodromic analogue is straightforward. Replacing d/dt —> 2irid/ dr, 
one obtains the non-autonomous system 

27^ = p„, 

2ni ~§^ = ~\ ^ a ' hup'ja ■ q)F- a F a , 
dF 

2ttz— ? = - J2 P(P • q)F a -pF p N a ,-p. (4.35) 
/3eA,Q-/3eA 

(Terms icluding G^s have been eliminated by the constraints.) These equations can be 
converted to the Lax equation 

wM + »|fl = [K,),M(,)]. (4.36) 

The monodromy of L(z) and M(z), too, takes the same form: 

L{z + l)=L{z), M(z + l)=M(z), 
L(z + r) = e 2niQ L(z)e- 2mQ , 

M{z + t) = e 2niQ {M{z) + 27TiL{z))e- 2niQ - 2mP, (4.37) 

where Q = X/u=i Qi*Hf*- The Lax equation implies that the monodromy data of the 
ordinary differential equation 

= L{z)Y{z) (4.38) 

on the torus E T is invariant as r varies. Y(z) now take values in the representation 
space V; the monodromy around a singular point or of a cycle of E T is represented by a 
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linear transformation on V. The ordinary differential equation has a regular singularity 
at z = only. The local monodromy around this singular point is a linear transformation 
T G GL(V). Similarly, the global monodromy along the a and (3 cycles give T a ,T/3 G 
GL(V). These linear transformations Tq, T a and are the monodromy data that are 
left invariant. 



5 Conclusion 

We have thus demonstrated that various models of the elliptic Calogero-Moser systems 
are accompanied with an isomonodromic partner. A technical clue is the choice of funda- 
mental functions x(u, z), y(u, z), etc. in the Lax pair L(z) and M(z). For L(z) and M(z) 
to give an isomonodromic Lax pair, these functions are required to satisfy a kind of "heat 
equation" besides the functional equations. We have illustrated the construction of the 
isomonodromic Lax pair for several typical cases — the Lax pair of the Ai_\ mode in the 
vector representation, the root type Lax pair for various untwisted and twisted models, 
and the Lax pair of the spin generalizations. 

The most interesting case in the context of Manin's equation is the root type Lax pair 
for the extended twisted BCe model (or, equivalently, the Inozemtsev system). The root 
type Lax pair based on short roots of the BCg root system consists of 2£ x 2£ matrices. 

The construction of a Lax pair, however, is merely the first step towards a full under- 
standing of Manin's equation and its possible generalizations. The next isse is to elucidate 
the meaning of the affine Weyl group symmetries, various special solutions, etc. in this 



framework. Recent works by Noumi and Yamada [2C], Deift, Its, Kapaev and Zhou [21] 
and Kitaev and Korotkin [[22[] are very suggestive in this respect. 

The spin generalization that we have discussed is a special case of a more general 
multi-spin system, i.e., the elliptic Calogero-Moser systems coupled to "Gaudin spins" 
sitting at the punctures of a punctured torus [§, |Hj. This is the Hitchin system on a 
punctured torus; we have considered the case with only one puncture located at z = 0. 
It is rather straightforward, though more complicated, to generalize our Lax pair to the 
multi-spin generalization. This gives a generalization, to other simple Lie groups, of the 



SU(2) isomonodromic system of Korotkin and Samtleben [23fl. The dynamical r-matrix 



in the work of Felder and Wieczerkowski |L7j plays a central role here. We shall report 
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this result elsewhere. 
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A Proof of Functional Identities and Heat Equation 
for Untwisted Models 

A.l Proof of (ggp 

Let f(u,v, z) denote the difference of both hand sides of fl2.8|): 

f(u, v, z) = x(u, z)y(v, z) - y(u, z)x{v, z) - x(u + v, z)(p(u) - p(v)). (A.l) 
This function turns out to have the following analytical properties: 

1. f(u,v,z) has the same quasi-periodicity as x{u, z) on the u plane, i.e., 

f(u + 1, v, z) = f(u, v, z), f(u + t, v, z) = e 2 "7(n, v, z). (A.2) 

2. f(u, v, z) is an entire function on the u plane. 

The first property is obvious from the quasi-periodicity of x(u, z) and the periodicity of 
p(u). Furthermore, poles of f(u, v, z) can appear only at the lattice points u = m + nr 
(m, ii G Z) on the u plane. Therefore, in order to verify the second property, we have only 
to show that f(u,v,z) is non-singular at these points. Actually, because of the quasi- 
periodicity, it is sufficient to consider the point u = only. As u — >• 0, the singular terms 
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x(u, z), y(u,z) and p(u) in f(u,v,z) behave as 

x(u,z) = - + 0(1), 
u 

y{u,z) = -A + 0(1) 



u 

p(u) = - 2 +0{u 2 ) (A.3) 



so that 



f{u,v,z) = (l + 0(l))y(v,z)-(-± + 0(l))x(v,z) 



U ■ . ■ \ u 2 



- (x(u, z) + y(u, z)u + 0{u 2 )) [ — - p(v) + 0(u 2 ) 



= 0(1). (A.4) 

We can thus verify the above two properties of f(u,v,z). 

Actually, any function with these two properties should vanish identically. This can 
be seen in several different ways. The shortest will be to resort to algebraic geometry 
of line bundles on the torus E T . A more elementary proof is to consider the quotient 
f(u,v,z)/x(u,z). This quotient is a doubly-periodic meromorphic function, and all pos- 
sible poles are located at the lattice points u = m + nr (m,n e Z), and at most of first 
order. In other words, f(u,v,z)/x(u,z) is a meromorphic function on the torus with the 
only possible pole at u = 0, but the order of pole cannot be greater than one. Such a 
function has to be a constant. On the other hand, because of the pole of x(u, z) at u — 0, 
f(u, v, z)/x(u, z) has a zero at u — 0. Therefore the constant should be equal to zero. 



A. 2 Proof of ( P^l) and flZTTU) ) 



(Of) can be readily derived from (|2.8[) by letting v — * — u. Let us consider (|2.10|). By 



(ID 



du 

Consequently, 



' x(u, z)x(—u, z) ) = — x(u, z)y(—u, z) + y(u, z)x(—u, z) = —p'(u). (A. 5) 



x(u, z)x(—u, z) = —p(u) + (independent of u). (A. 6) 
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Since x(u, z) = —x(z,u) = —x(—u, —z), the left hand side of the last relation is in fact 
an anti-symmetric function of u and z. Therefore, 



x(u, z)x(—u, z) = p(z) — p{u) + const. (A. 7) 

Now consider the limit as u — > z. Both x(u, z)x(—u, z) and p(z) — p{u) tend to zero in 
this limit. Thus the constant on the right hand side has to be zero. 



A.3 Proof of (1Z7T81) 



Let us rewrite the both hand sides of fl2.18|) into a more accessible form. Differentiating 
x(u, z) by r gives 

= x^z^ilogO^z -u) + log^(O) - \o g e 1 (z) - lo g e 1 (u)). (A.8) 
By the heat equation ( 2.24|) of the Jacobi theta function, 



d n , , 9"(u) d [6\(u)\ fe'Au) . 
4ttz— 6i(u) = = IT TrH + 7TTT = P ( u ) + P( u ) • ( A - 9 ) 



dr 1V ' d x {u) du \9 x (u) J \0i(u) 
Letting u — > and recalling the singular behavior of p(u) at u = 0, we obtain 

47ri- log^(0) = hm(p'H + ?(«))= (A.10) 

Plugging these formulae into the above expression of dx(u, z)jr gives 

4m^A = x(u,z)f(u,z), (A.ll) 

where 

f(u, z) = p'(z -u) + p(z - uf + J^i - p'(^) - p(^) 2 - p'(u) - p(uf. (A.12) 

On the other hand, we have 

dx(u, z) d 
dudz 

where 



~dz Z ^ P ^ + p{yZ ~ = ~ X<kU ' z ^ 9 ^ z ^ (A.13) 



g(u, z) = (p(z -u)- p(z))(p(u) + p(z - u)) + p\z - u). (A.14) 

The goal is to verify that f(u, z) = 2g(u, z). It is sufficient to prove the following two 
properties of f(u, z) — 2g(u, z), because such a function has to be identically zero. 
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1. f(u, z) — 2g(u, z) is a doubly-periodic function on the u plane with primitive periods 



2. f(u, z) — 2g(u, z) is an entire function, and has a zero at u — 0. 

The first property is obvious if one notices the following quasi-periodicity of f(u,z) and 
g(u,z): 

f(u + 1, z) = f(u, z), f(u + r,z) = f(u, z) + im(p(u) + p(z - u)), 

g(u + l,z) = g(u,z), g(u + t, z) = g(u,z) + 2iri(p(u) + p(z -u)). (A.15) 

Let us check the second property. Possible poles of f(u, z) and g(u, z) are located at the 
two points u = and u = z of the fundamental domain of the period lattice Z + rZ. 
Again recalling the singular behavior of p(u) at u — 0, one can confirm by straightforward 
calculations that 



Thus f(u, z) — 2g(u, z) turns out to be non-singular and have a zero at u = 0. Similarly, 
one can see that f(u, z) — 2g(u, z) is non-singular at u = z. 



B Verification of Lax Pair for Extended Twisted BCi 



To prove the Lax equation, it is sufficient to derive the following three equations: 



1 and t. 



f(u,z) = 0(u), g(u,z) = 0(u) (u 



0). 



(A.16) 



Model 



dX a (z) 
dt 

dp ■ /i 



[P,X a (z)} (a = 1,2, 



[X x {z) + X 2 {z) + X 3 (z), Y 1 (z) + Y 2 (z) + Y 3 (z)U, 
[X x {z) + X 2 (z) + X 3 (z), D + Y l {z) + Y 2 (z) + Y 3 (z)} 



3), 



(B.l) 



(B.2) 



dt 





(/i^).(B.3) 



p and v run over the set A s of short roots. 

The proof of ( |B.1| ) is quite easy. Let us consider the case of a = 1. The t-derivative of 
Xi(z) can be written 




(B.4) 
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Using the commutation relation [P, E(a)] = a ■ pE(a), one can readily see that the right 
hand side is equal to [P, Xi(z)]. The other two in ( |B.1| ) can be similarly derived. 
The rest of this appendix is devoted to the other two equations ( P-2| ) and (|B.3| ) . 



B.l Proof of (\B72j ) 

We calculate the diagonal elements 

[X a (z), Y b {z)}„ = [x a ^{z)Y b ^{z) - Y b ^{z)X a ^{z)) (B.5) 
of the nine commutators one-by-one. 
B.l.l Vanishing terms 

Some part of the matrix elements of X a (z) and Y b {z) turn out to vanish by the nature of 
the BCi root system: 

Xi,^(z) = Y x „{z) = 0, (B.6) 
X 2 ,^(z)=Y 2 ^(z)=0 (ji^-u), (B.7) 
X 3 ^(z)=Y 3 ^(z)=0 (ji^-v). (B.8) 

The first relation is due to the fact that fi — (— /i) = 2fi can never be a middle root. 
The second and third relations are obvious if one notices that — v is a long root (or, 
equivalently, twice a short root) if and only if [i = —v. 
In particular, 

\X x {z),Y 2 {z)\ m = y 3 (*)W = [M*),Yi(z)U = IM^M^U = 0. (B.9) 

B.l. 2 Calculation of [X^z), 

By definition, 

[x 1 (z),Y 1 (z)] (tll = -gfa Yl ~ ^ ' 9 > z }v(( v ~ d ' q > z> > 

-y{{v z)x{(v -fj,)-q,z)). (B.10) 
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We rewrite this sum to a sum over the middle root a = // — v. Since the middle roots a 
of this form are characterized by the condition that a ■ \i = 1, the right hand side can be 
rewritten 

agA m ,a-/i=l 

Actually, the possible values of a ■ [A are limited to and ±1 only. Therefore this sum is 
equal to 
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a ■ fj,(x(a ■ q, z)y{—a ■ q, z) — y{a ■ q, z)x(— a ■ q, z)\ . 



(The factor 1/2 compensates the contributions from a ■ fi = 1 and a ■ fi = —1.) Noting 
that a ■ fi = {p ■ fi, a ■ q} , we can express [X^z), ^(-z)] as a Poisson bracket of the form 

[X 1 (z),Y 1 (z)]„ i = {p-li,V n }, (B.ll) 

where 

Vu = ^f Yl < a ■ 9. z X~ a • ( R12 ) 

B.1.3 Contributions of other commutators 

By ([B.7|) and (|B.8|), the diagonal elements of the other commutators are a sum of just 



two terms: 

[X a {z), Yf,(z)]^ = Xa^-fjXb,-^ ~ Yb^-^Xa^^. (B.13) 

Let us consider the case of a = 2 and b = 2 in some detail. By definition, 

[X 2 (z),Y 2 ( z )U 

= ~ (gixx&H -q,z) + gi2X {2) (2/i -q,z)) (fl^y (-2// -q,z) + g [2 y {2) (-2/i • q, z) J 
+ [9iiy(^ ■ q,z) + gi 2 y (2) (2fi ■ q,z)j (g a x(-2n ■ q, z) + g l2 x (2) (-2fi ■ q,z)j. 

Since a = 2/x is a long root, and long roots with non- vanishing inner product with fi are 
2/x and — 2/x only, the right hand side can be rewritten 



~\ a ' A*(fti«(a • ?>z) +5'/2a; {2) (« • (gny{-oc • q,z) + g i2 y {2 \-a ■ q, ; 

+ i a ' A 4 (fltaj/fa -q,z) + 9i2y m {oc ■ q, z)j (gnx(-a -q,z) + gi 2 x {2) (-a ■ q, z)j . 



aGA; 
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(The factor 1/4 compensates the contributions from a ■ [l — 2 and a ■ /i = —2.) We can 
again cast this into a Poisson bracket: 

[X 2 (z),Y 2 (z)]^ = {p-^V 22 }, (B.14) 

where 

Vti = 4 S ' 9' 2 ) + 5" 2x(2) ( a ' 9' z )) -Q,z) + 9i2X {2) (-a ■ q, zfj . 

aeAi 

(B.15) 

Similarly, one can obtain 

[X 2 (z), Y 3 {z)U = W, V23}, [Xs(z), Y 2 (z)U = {p ■ M) F32}, 
[X3(^),r3(^)U = {p-/i,^33}, (B.16) 

where 

V23 = ^ ^2 (gn x ( 2a -q,z) + gi 2 x (2) (2a ■ q, z)j (g sl x(-a ■ q, 2z) + g s2 x {1/2) (-a • q, 2zfj 

o£A s 

V32 = - ^ ^ S ix(a • g, 2z) + g s2 x {1/2) (a ■ q, 2z) s j (g a x(-2a -q,z) + g i2 x {2) (-2a ■ q, z)j 

«£A S 

V33 = ^2 (gsix(a ■ q, 2z) + g s2 x (1/2) (a ■ q, 2z)j [g s ix(-a ■ q, 2z) + g s2 x {1,2) (-a ■ q, 2z 



a€A s 



(B.17) 



Collecting the results of these calculations, we find that the right hand side of ( |B.2| ) 
takes the form of the Poisson bracket {p ■ fi, V}, where 

V = V n + V 22 + V 23 + V 32 + V 33 . (B.18) 
B.1.4 Writing V in terms of p functions 

The final step is to rewrite V in terms of the Weierstrass p functions. For Vu, this can 
be done by use of (|2.10|) . The other parts are due to the following functional identities: 



x^ 2 \u,z)x^ 2 \-u,z) = -p^ 2 \u) + p (1 / 2 \^), (B.19) 
x {2) (u,z)x {2) (-u,z) = -p {2) (u) + p {2) (2z) } (B.20) 
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x(u,2z)x {1/2) (-u 1 2z) + x (1/2) (u,2z)x(-u,2z) = -2p(u) + const., (B.21) 

x(u, 2z)x(—2u, z) + x(2u, z)x(—u, 2z) = —p(u) + const., (B.22) 

x(u, 2z)x {2) (-2u, z) + x (2) (2u, z)x(-u, 2z) = -p{u) + const., (B.23) 

x {l/2 \u,2z)x(-2u,z) + x(2u,z)x {1/2) (-u,2z) = -p {1/2) (u) + const., (B.24) 

x (1/2) (u, 2z)x {2) (-2m, z) + x {2) (2m, z)x {1/2) {-u, 2z)} = -p(u) + const., (B.25) 

x(u, z)x^ 2 '(—u, z) + x^ 2 '(u, z)x(—u, z) = —2p( 2 \u) + const. (B.26) 

The first two are substantially the same as ( [2.10| ) except that the variables and the 
primitive periods are rescaled. "const." in the other identities stand for terms that are 
independent of u, thereby negligible in the Poisson bracket with p- «; remember that they 
are not absolute constants, but functions of z and r. We shall prove these identities in 
Appendix C. Using these functional identities, one can see that V is equal to the potential 
part of the Hamiltonian 7i, up to non-dynamical terms independent of p and q. 

To summarize, we have shown that the sum of the (w, w) elements of the nine com- 
mutators coincides with the Poisson bracket {p ■ u, V}, which is equal to dp ■ fi/dt by the 
equations of motion of the model. 

B.2 Proof of (IB73I ) 

The proof can be separated into the cases where v = — u and v ^ ±u. 
B.2.1 u = -fi 

The vanishing of the (w, — w) elements of the commutators other than [X a (z),D] (a = 
1, 2, 3) and [Xi(z), D] is immediate from (|B.7| ) and (B.8| ). [X a (z), D]^ _^ vanishes because 
of the symmetry _D_ M = D^. As for [Xi(2;), Y^z)]^^, we have 

[Ai(z),ri(2:)U- M = -g 2 m x (^-v) ■ q, z)y((v + -q,z) 

+9 2 m J2 y((t*-")-q,zM(v + tJ>)-q,z)- (b.27) 

veA 3 \{±n} 

Bu substituting v — > —u, the second sum on the right hand turns out to be identical to 
the first sum. The two sums thus cancel with each other. 
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B.2.2 v^±fi 



The following can be readily seen by using (|B.7|) and ( 



[X 2 {z),D] l „,= [X 3 {z),D] lu , = Q, 

[X 2 (z), Y 2 (z% v = [X 2 (z), Y 3 (z)U = [X 3 (z), Y 3 (z)l v = 0. (B.28) 
The (/i, u) elements of other commutators can be calculated as follows: 



[Xi (*),£>] 



[Xi{z),Y>{z)\ 



[X*{z),Y x {z)\ 



—Xi ; ^(z)(D fl — D h 
g m x({p -u)-q, z) 



x \9aip(fi ■ q) + g s2 p (1/2 \fx ■ q) + #np(2/i • q) + gi 2 p (2) (2fx ■ q) 
-9sip(y ■ <?) - g S 2p {1/2) (v ■ q) - gnp(2u ■ q) - g l2 p {2) (2v ■ q) 



+ Y p( a -<i)- Y p( a -i))- 

AGA m ,o-^t=l ct(zA m ,ofU=l 



(B.29) 



Y (x^ x (z)Y 1M (z) - Y ltilX (z)X lM (z) 
\ga s 

\eA s \{iJ,,u} 
-y((/i - A) • q, z)x{{\ -v)-q, z)\ . 



(B.30) 



X l 

—v,v (z)Xi, -\x,v\Z) 

-g m x((ii + v) ■ q, z) (gny(-2v -q,z) + g i2 y {2 \-2u ■ q, z)j 

+ (fl'/iy(2/i • q,z) + g l2 y {2) (2/x • q, z))g m x(-(fj, + v) -q,z). 

X! (z)Ys- U)l ,(z) — Y3 t n- ll (z)Xi- lJt)V (z) 

-g m x((ii + u)-q, z) \g aX y(-v ■ q, 2z) + g S 2y {1/2) '(-f • q, 2z) 

+ (gaivip ■ ?, 2z ) + g S 2y il/2) (n ■ q, 2zyjg m x(-{n + v) ■ q, z). 

X 2 ^- ll {z)Yx,- t i,v( z ) — Yi tfM - u (z)X 2 - u ,u( z ) 
-(g n x(2fj l ■ q, z) +gi 2 x {2) {2^ ■ q, z)^g m y(-(n + v) ■ q, z) 



(B.31) 



(B.32) 
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+g m y((fi + v)-q, z) \gnx{-2v -q,z)+ gi 2 x {2 \-2u ■ q, z)j . 

(B.33) 

[X 3 (z),Y x (z)U = X 3 ^ v (z)Y 1 ^{z)-Y 1 ^{z)X 3 ^ >u (z) 

= -2^g sl x(fi ■ q, 2z) + g s2 x {1/2) (p, ■ q, 2z)^g m y(-(fi + u)-q, z) 
+2g m y((/j l + v)-q, z) (g sl x{-v ■ q, 2z) + g s2 x {ll2) (-v ■ q, 2z)\ . 

(B.34) 

We now sum up all these quantities, regroup terms into those multiplied by the same 
monomial of coupling constants, and show the cancellation in each partial sum. There 
are six monomials of coupling constants that can occur — g^, g m gn, g m gi2, 9m9s\ and 

9m9s2- 

Let us consider the terms multiplied by g^. This is a sum of the following two quan- 
tities: 

I = x((n-v)-q,z)( Y P( a -<l)- Y Pia-V)) 

ct£A m ,a-/^=l a&A m ,a-v=l 

II = - Y ( x ((^ - A ) ■ ?>%(( A - v) ■ q,z) 
-y((j* - A) ■ g, z)x((X -u)-q, z)j . 
By the functional identity (|2.8| ), we can rewrite II into a sum over middle roots: 

ii = - y a; ((/ i - z/ )-^ 2; )(p((/ i - A )-9)-p(( z/ - A )-?)) 

\£A s \{fi,v} 

= -x((fx - v) ■ q, z) ( Y P( a -q)~ Y 

Here the sum over A has been converted to a sum over a by putting a = fi — A and 
a = v — A in the two p function in the first line. Note that /x, v and A are all orthogonal 
to each other. We thus find that 1 + 11 = 0. 

For the other partial sums, we use the following functional identities, which we shall 
prove in Appendix C: 

x(2u, z)y(—u — v, z) — y(2u, z)x(—u — v, z) + x(u + v, z)y(—2v, z) 

—y(u + v, z)x(—2v, z) — x{u — v, z)(p(2u) — p(2v)) = 0, (B.35) 
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x^ 2 \2u ) z)y(—u — v, z) — y( 2 \2u } z)x(—u — v , z) + x{u + v, z)y ( - 2 \—2v, z) 

-y(u + v,z)x {2 \-2v,z) - x(u-v,z)(p {2) (2u) - p [2) (2v)) = 0, (B.36) 

2x(u, 2z)y{—u — v , z) — y(u, 2z)x{—u — v, z) + x{u + v, z)y(—v, 2z) 

—2y{u + v, z)x(—v, 2z) — x(u — v, z)(p(u) — p(v)) = 0, (B.37) 

2x (1/2) {u, 2z)y(-u -v,z)- y (1/2) {u, 2z)x{-u -v,z)+ x{u + v, z)y {1/2) (-v, 2z) 

-2y(u + v, z)x {1/2) (-v, 2z) - x(u - v, z)(p (1/2 \u) - p (1/2) (v)) = 0. (B.38) 

By these functional identities, we can confirm that all the partial sums regrouped by 
g m gn, g m 9i2, g m g s i and g m g S 2, respectively, cancel out. 



C Proof of Functional Identities for Twisted Models 

We here prove the functional identities that we have encountered in Appendix B. Although 
the proof is optimized to our choice of x(u, z), x^'^'iu, z) and x^ 2 '(u, z), the same method 
can in principle apply to other solutions of the functional equations, such as the functions 



used by D'Hoker and Phong [15] and Bordner and Sasaki [14 . 



C.l Analytical properties of x^^ 2 \u, z) and x( 2 \u, z) 

The proof of the identities including x^^ 2 \u, z) and x^(u, z), like the proof in Appendix 
A, is based on the analytical properties of those functions. 

• x^ 1 ' 2 ^ (u, z) has the following analytical properties: 

1. x^^ 2 \u, z) is a meromorphic function of u and z. The poles on the u plane and 
the z plane are located at the lattice points u = m/2 + nr and z = m + 2nr 
(m, n G Z). 

2. x^^ 2 \u, z) has the following quasi-periodicity: 

x^ 2 \u+-,z) =x^ 2 \u,z), x^ 2 \u + r,z) = e 2mz x^ 2 \u,z), 

(u,z + l) =x {1/2) (u,z) } x il/2) (u,z + 2r) = e^ iz x {1/2) (u,z).(C.l) 



X 



(1/2) 
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3. At the origin of the u and z planes, this function exhibits the following singular 
behavior: 

x (1/2 \u, z) = -- 2p{z | 2r) + OOu) (u -> 0), 
u 

x (1/2) {u,z) = -- + 2p{2u\2r)+0{z) (z -> 0). (C.2) 
x( 2 '(u,z) has the following analytical properties: 

1. x( 2 '(u, z) is a meromorphic function of u and z. The poles on the u plane and 
the z plane are located at the lattice points u = 2m + nr and z = m + nr/2 
(m, n£Z). 

2. x( 2 )(w, z) has the following quasi-periodicity: 

x (2) (u + 2, z) = x (2) (u, z) , a; (2) (u + r, *) = e 27r "x (2) (u, z) , 
x {2) (u,z + l) = x {2) (u,z), x {2 \u,z + ^)=e mu x^ 2 \u,z). (C.3) 

3. At the origin of the u and z planes, this function exhibits the following singular 
behavior: 

x^(u,z) = -- 1 -p(z\^) + 0(u) («-0), 
it 2 2 

x^(u 1 z) = -j- + 1 - p q\^)+0(z) (z-0). (C.4) 



C.2 Proof of (Em) - ( Egg ) 



These four identities can be treated in much the same way. Let us illustrate the proof for 
( |B.35p only. Since the line of the proof is almost the same as the proof of ( |2.£| ), we show 
an outline of the proof and leave the details to the reader. 
Let f(u,v,z) denote the left hand side of ( p.35|) : 

f(u,v,z) = x(2u, z)y(— u — v, z) — y{2u, z)x(— u — v, z) + x(u + v, z)y(— 2v, z) 

—y{u + v, z)x(—2v, z) — x{u — v, z)(p(2u) — p(2v)). (C.5) 

Our task is to show the following analytic properties of f(u,v, z), which imply that this 
function is identically zero: 

45 



1. f(u,v,z) has the quasi-periodicity as follows: 

f(u + 1, v, z) = f(u, v, z), f(u + r, v, z) = e 2mz f(u, v, z). (C.6) 

2. f(u, v, z) is an entire function on the u plane. 

The first property is immediate from the quasi-periodicity of x(u, z), etc. Furthermore, 
it is obvious from the definition that all possible poles of f(u,v,z) on the u plane are 
limited to the lattice points u = m/2 + nr/2 and u = —v + m + nr (m, n G Z). In view 
of the quasi-periodicity, therefore, we have only to verify that f(u,v, z) is non-singular at 
u = 0, 1/2, r/2, 1/2 + r/2, and -v. 

The absence of poles at u = 0, 1/2 and — v can be verified by straightforward cal- 
culations on the basis of the singular behavior of x(u,z), x^^ 2 ^(u,z) and x^(u, z) as 
u — > 0. 

In order to examine the points u = r/2 and u = 1/2 + r/2, one has to examine the 
singular behavior of x(2u, z) and y(2u, z) as u — > r/2, 1/2 + r/2. This can be worked out 
by combining the quasi-periodicity of x(u, z) and y(u, z) and their singular behavior as 
u — > 0'. 



1. As u -> r/2, 



x(2u, *) = e 27r ^x(2M - r, z) = e 27riz ( — *— + 0(1) ) , 



y(2u,z) = 2 2mz y(2u-T,z)=e 2mz (- /n 1 .„ + 0(1) ) ■ (C.7) 



2m -r 
1 

(2^r) 2 



2. As « -> 1/2 + r/2, 



x(2u, z) = e 2 ^ z x(2 M - 1 - r, 2) = e 2 ^ 2 ( + 0(1) ) , 



y(2u,z) = e 2mz y(2u — 1 — r, z) — e 2mz I — -^ + 0(1)). (C.i 



2u - 1 - r 
1 

(2m - 1 - r) 5 

Using these observations, one can confirm the absence of poles of f(u,v,z) at u = r/2 
and 1/2 + r/2 by direct calculations. 

We can thus verify that f(u, v, z) is indeed an entire function on the u plane. 
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C.3 Proof of (07ZT\ ) - ( ggg) 

Rather than directly proving these identities, let us prove them in a differentiated form. 



For illustration, we consider the first identity ( ]B.21|) . Differentiating this identity by u 
gives 

x(u, 2z)y ( - 1/2) (-u, 2z) - y{u, 2z)z {1/2) (-«, 2z) + x {1/2) (u, 2z)y{-u, 2z) 
-y {1/2) (u, 2z)x(-u, 2z) = 2p\u). (C.9) 

One can prove it directly, repeating the complex analytic reasoning that we have presented 
in other cases. An alternative way is to take the limit, as v — > u, of the functional identity 

x(2u, 2z)y {l/2 \-u - v, 2z) - y(2u, 2z)x {1/2 \-u - v, 2z) + x {l/2) (u + v, 2z)y(-2v, 2z) 
-y {1/2 \u + v, 2z)x{-2v, 2z) - x{u - v, z){p{2u) - p{2v)) = 0. (CIO) 

(This yields the above identity upon substituting u —>■ u/2 and v — ► v/2.) This functional 
identity can be derived by the same method as the proof of (p.35| ) - ( p.38|) . 



Similarly, the third and fifth of QB.21 ) - (p.26| ) are obtained from the following func- 
tional identities: 

2x(u, 2z)y^ 2 \—u — v, z) — y(u, 2z)x^ 2 \—u — v, z) + x^ 2 \u + v, z)y(—2v, 2z) 
-2y {2 \u + v, z)x(-2v, 2z) - x(u - v, z){p{u) - p{v)) = 0, (C.ll) 
2x {1/2) (u, 2z)y (2) (-u -v,z)- y (1/2 \u, 2z)x {2 \-u -v,z) + x {2 \u + v, z)y (1/2) (-2v, 2z) 
-2y ( - 2) (u + v,z)x ( - l/2) (-2v,z) - x(u-v,z)(p(u) - p{v)) = 0. (C.12) 

The second, forth and sixth of flBjjp - (|BT26|) can be similarly derived from the last 
three of ( |B.35|) - ( |B.38|) . This completes the proof of the functional identities. 

We conclude this appendix with a comment on the "const." terms of these identities. 
In principle, these terms can be determined by examining the identities at a special point 
of the u plane. Let us consider, e.g., ( [B.21|) . At u = z, the first term on the left hand side 
vanishes. Evaluating the other terms at this point, therefore, one finds that 

const. = 2p(z) - x {1/2) (z, 2z)x{-z, 2z). (C.13) 

The same formula can be reproduced by substituting u = —z. One can similarly derive 
an explicit expression for the other identities. 
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